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We derive analytically the asymptotic behavior of the Casimir interaction between a sphere and 
a plate when the distance between them, d, is much smaller than the radius of the sphere, R. The 
leading order and next-to-leading order terms are derived from the exact formula for the Casimir 
interaction energy. They are found to depend nontrivially on the dielectric functions of the objects. 
, As expected, the leading order term coincides with that derived using the proximity force approx- 

imation. The result on the next-to-leading order term complements that found by Bimonte, Emig 
' and Kardar [Appl. Phys. Lett. 100, 074110 (2012)] using derivative expansion. Numerical results 

are presented when the dielectric functions are given by the plasma model or the Drude model, with 
5_( . the plasma frequency (for plasma and Drude models) and relaxation frequency (for Drude model) 

given respectively by 9eV and 0.035eV, the conventional values used for gold metal. It is found that 
if plasma model is used instead of Drude model, the error in the sum of the first two leading terms 
is at most 2%, while the error in Oi, the ratio of the next-to-leading order term divided by d/R to 
the leading order term, can go up to 4.5%. 
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I. INTRODUCTION 



Casimir effect is a quantum effect that cannot be ignored in the realm of nanotechnology. It can cause malfunctions 
of nano devices due to stiction [H-Q. In the last decade, intensive research have been carried out to determine the 
exact analytic formula for the Casimir effect between two nonplanar objects and its effective numerical computations 
(see, for example, the references cited in Q). Prior to this, one can only rely on the proximity force approximation 
(PFA) to compute an approximation for the Casimir interaction, and there is no way to determine the magnitude of 
. the error in such an approximation. 

In the case of the sphere-plate setup, the most popular configuration used in Casimir experiments, there is only 
' one curvature parameter given by the radius of the sphere, R. Hence it is expected that as d, the distance from the 
' sphere to the plate, is much smaller than i?, the Casimir interaction energy has an asymptotic expansion of the form 



Ec., = E^l^(l + ^e,,E + ...]. (1) 



JH ' where E^^^ is the proximity force approximation to the Casimir interaction energy. It follows that for the Casimir 
- -' force Fcas and force gradient dFca.s/dd, one also has expansions of the form 



-dT-^d-V^R^' 



(2) 



A few years ago, experiments have been set up to measure 9i using a micromachined torsional oscillator (5|. This gives 
a more ernest reason for the theoretical computation of the next-to-leading order terms of the Casimir interaction. 
One of the breakthroughs in Casimir research brought by the achievement in explicit functional representation of the 
Casimir interaction is that it becomes possible to compute analytically the next-to-leading terms, as has been shown in 
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for the cylinder- plat e configuration, in [9l-[ll| for the sphere-plate configuration, in for the cylinder- cylinder 
configuration, and in [ij] for the sphere-sphere configuration. However, except for Q, all the other works only deal 
with ideal or non-physical boundary conditions, i.e., Dirichlet, Neumann, perfectly conducting, infinitely permeable 
or Robin boundary conditions. So far no work has discussed the exact analytical computation of the next-to-leading 
order term in the Casimir interaction between a sphere and a plate when both of these objects are made of real 
materials, and this is the goal of the current work to deal with this problem. 

It should be mentioned that there has been an attempt to compute the material dependent next-to-leading order 
term in the Casimir interaction between a sphere and a plate carried out by Bimonte, Emig and Kadar which 
used the method of derivative expansion postulated in fisi. which in turn is inspired by the work [l^. However, it 
is still desirable to check the validity of the postulate in |l4 [T5| by computing the next-to-leading order terms from 
the exact formula for the Casimir interaction. Therefore, the results of our current work complement those obtained 




II. THE CASIMIR INTERACTION ENERGY 



In this article, we recall the formula for the Casimir interaction between a sphere and a plate. Assume that the 
sphere has relative permittivity e^.i, and the plate has relative permittivity £^.2- When the thicknesses of the sphere 
and the plate are larger than their respective skin-depths, we can model this configuration by a ball and a semi-infinite 
space. Let d be the distance from the sphere to the plate, and let L = d + R, where R is the radius of the ball. 

As shown in [itI [isj , the electromagnetic Casimir interaction energy of this sphere-plate configuration is given by 



Ecs = ^ / deTrln(l-M(zO), 



where the trace Tr is 



(3) 



Tr 



m=0 i=max{l,|m|} 



with tr being the trace over 2x2 matrices. The matrix elements of M are given by 

Hm,l'm' =0m,„,' M . — 7777^ 77 Ji / d9 Smh 9 e'^"^ """""^ ^ 



+ + 1) (/ + m)!(/' + m)! 



sinh0 Pi"" {cosh 9) 
-Pr {cosh 9) 



sinh 



ihfi^^"(-^^^)' 
sinh 6IP,™' (cosh 61) 



^ sinh 9 PlI''' {cosh 9) 



-P7^ (cosh0) 



m 



sinh0 



PIT (cosh 6*) 



sinh0 

sinh 9 pit' '{cosh 9) , 



where 



is a diagonal matrix with elements 



^™(^0 



/TiTE n 

^ -I t™ 



Ii+i{iiR) (niKi?) -1- niKRI[^ 


1 {niKR)j - Ii^i{niKR) (|/;+ 1 {kR) 4 


- KRI[^^{nR)^ 


Kij^i{kR) 1 (ni^i?) + uikRI'^^ 


i{niKR)j — (niAti?) (^^Ki_^_i{kR) - 


V hRK[^,_{kR)^ 



Ki^i{kR) (^^Il_^^l{nlKR) + niKRI[^^{niKR)^ - £^,1/,+ 1 (ni^i?) (^\Ki^i{kR) + kRK[^^{k,R)^ 



and 



T™ ^0 
T™ 



(4) 
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is a diagonal matrix with elements 



-La — 



+ sinh^ 9 — cosh ( 



i2 I o^T->l->'^ 



sinh 6 + cosh L 



sinh^ — er2 cosh 



sinh^ 6* + e,. 2 cosh ( 



Here 



-i - ^ --19 



C 

and Pl^{x) are the associated Legendre functions given by 



when m > 0, and 



Direct numerical computations of the Casimir interaction energy from the formula have been performed in a 
few works, for example, in (l9l . [20| . In numerical computations, the infinite matrix M has to be truncated to a 
matrix of finite size. A drawback of this direct numerical computation is that when d/ R gets smaller, one has to 
use a truncated matrix of larger size for accuracy, and this is subjected to the capacity of the computer. Currently, 
numerical computations are limited to d/R > 0.05. However, in experiments, we usually have d/R ~ 0.01. Hence, 
analytical computation of the Casimir interaction energy becomes desirable. 

III. SMALL SEPARATION ASYMPTOTIC EXPANSION 

In this section, we want to derive analytically the small separation asymptotic expansion of the Casimir interaction 
energy, Casimir force and the force gradient up to the next-to-leading order term. 

One of the technical issues in the analytical computation of the Casimir interaction energy ^ is the appearance of 
the associated Legendre functions P"'(a;). First notice that because of the relation ^ and 

I -2k — \ /, „ N / / - 2fc 



sinh6' (1 \ I sinhe W 1 , . . 



sinhfl ' ^ sinh0 

the matrix element M.im,im Q is equal to that when m is changed to — m. Hence, it is sufficient to consider nonnegative 
TO. In this case, one can show that 

P,'"(coshg) =f-iynim il + ^)l I ^(^(-cogi^5/_^sinh6lcos^)'cosTO^ 

^(_^),n^,n([+Z^y 1 n-2k)e n d^cos2'-2fc(psin2'=^e2™^. 
TT ^-^ k\{l — k)\ J 

Differentiating with respect to gives 

sinh (cosh 0) ={-iri^ + m)l ' I - 2k n^^^^ H rf^ cos^'-^'^ ^ sin^'^ ^e^™^. 

TT ^-^ k\(l — k)\ / TT 
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Therefore, 



-F;" (cosh 61) 



/smh6iP™'(cosh6l) 
I smh^ 

I 777 

\ ——p,'" (cosh 61) sinh6'P,""(cosh6') 



sinh 9 



TT ^^k\(l-k)\ 



sinh 



l~2k, 



Making a change of variables 

and expanding the logarithm in ([3]), we have 

he 



Cas 



2ttR 



OO ^ „oo OO / S OO \ / S 

Et^t/ '^-e n E 

s=0 "^0 m=0 yi=0/i=maxl,|m|/ \i=a 



where 



= — V(2^^ + l)(2Z,+i + 1)(Z, - m)l{h+i - m)\{h + my.{k+i + my. i 



-T™ 



U h 



EE 



fe=0 fe'=0 



fc!(Z, -A:)!fc'!(/,+i-A:')!7o 



^6* sinh 6'e"^'^(^+^) """"^ e+{u+u+t~2k~2k') 



I k~2k 



m L — 2k 



2k' 



~T™ 



V^.+i(/»+i + l) 
m 



sinh 6* + 1) 

- 2fc' 

Vsinh6'v/^»+i(/»+i + l) V^»+i(^i+i + 1) ) 



d^cos^'^-^'^^sin^'^^pe^™^ 



VcOS^'^+i-^fc' ^'sij,2fc' ^'g2W^ 



(7) 



with 



4+1 (w) (^i/;^+i(niw) + niw/;^_^i(niw)^ 


1 -a^/,,+ i(niw) 


(i/,.+ i(^)+<+iM) 




- "t^+i ("1^) (i^i.+i (^) + (^)) 



Here * — TE or TM, and a™ = 1, a™ = e^.i- The minus signs on Tj"^ and Tg^^ in ([7]) come from the two matrices 



1 

-1 



in ©. Let 



In the following, we make a shift of parameters 

I— ^ / + /i, 

where 



sinh = —. 

UJ 
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When e ^ 1, the main contributions to the Casimir interaction energy come from the terms with 



7 1 7 1 1 1 fl 



In the small e expansion below, we will count the order oi 1,^,111,10 and as 1 /e, l/\/e, 1 1 /e and e respectively. 
Making a change of variables 



UJ = 



we have 

where Iq — and Is+i = by default. The integration over 6 is from — to 00 which can be approximated by an 
integration from —00 to cxi, since 9 is of order e and 9q is of order 1. 

Now we perform the small e expansion of Writing cosip as exp (— In sec (p) and using the fact that 



we have the following small e expansion: 



(^2 <^4 

lnsec^=- + -, 



2k 



d^cos^C+'-'-^'^^sin^'^^e^™^ « / ^ d^^^fe ( 1 - ^ j exp (^-{l + k- k)ip^ - i±I^^^ ] e^™^ 

1 r , 2k f, W\ ( l + k~k ^ l + h-k . 2imip\ 

ztmip 



J d(^¥'"'(l + A,2)exp(6,,i+6,,2)exp(^-^2 ^ , ^ 



(9) 



In the second line, we have performed a rescaling if 1— tp/Vl so that the main contribution to the integration over (p 
comes from ip that are ~ 1. Here and in the following, for any X , Xi^i and Xi^2 are, respectively, terms of order ^/e 
and e. When these terms do not depend on j, i would be omitted. Changing U to U^i and fc to fc' in (|9]), we obtain 
a similar expansion: 

j\ d(p'cos2('+''+^)-2'=' ^'sin^'^' ^'e^""^' V^'^' (1 + C,^) exp + P,,^) exp (^-^'^ + ^) . 

Next, we can use Stirling's formula 



In n! = \ n + - ) In ?i — n + - In 27r + — — 

2 12n 



to obtain an expansion 



1 ^J{1 + h- m)l{l + k+i - m)\{l + k + ui)\{l + k+i + m)! /m^ 

(i + k-k)\{i + k+,-ky. " {— + 

On the other hand, we have 

^x/(2; + 2k + 1){21 + 2h+i + 1) w (1 + + ^,,2) . 



For the terms involving 9, expanding in small e gives 

,tln fl„ -I 1 ^ 

VI 



02 \ T 

sinh(6l + 6I0) « sinh6io ( 1 + ^coth^o + ~ , (1 + + £^,2) ; 
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exp(^-2a;(l + s) cosh(6l + 60) + {21 + + k+i -2k- 2k'){e + 60)^ 
> exp((2Z + + k+i -2k- 2k')eo) 

-2w(l + £)sinh6'o f coth6lo + 9 + — cotheio + y + 24 '^^^'^^o ) + + k + k+i -2k- 2k')9 



1 + r 



k+k'- 



exp 



21 16^ 2el n , ,n ^ 

h [h + H+l)9 + J-i,l + J-i,2 I ; 

T T T 



l+li-2k 



m 



+ + + smh{e + eoW{i + ii){i + ii + i) 

m I + li — 2k 



\ 



I + h+i - 2k' r 



I 



1 + A,2 Mx 

Mx 1 + A,2 



V'(/ + ?i+i)(? + ^,+i + l) 
m 



sinh(0 + 0o)V(' + ^i+i)(^ + ^i+i + l) 
/ + /i+i - 2k' 



\ 



1+M,2 My 



Here 



Ml = 



mVl — 



It 



is of order We do not need the term that is of order e for the off-diagonal terms of these matrices as they won't 
contribute to the next-to-leading order term of the Casimir interaction energy. Finally, the small e expansions of 
^0+00 same as the small expansions: 



rri* 

^e+00 - 



^Jn^+s^Dh^{e + 60) - al cosh(6' + 9o) 
^Jnl+smh^{e + 60) + 0*2 cosh.{e + 9o) 



where aj^ = 1, aj"^ = er,2, sgn'^'^ = 0, sgn^^^ = 1, 

V£r,2(l-r2)+r2-l 



rriTE 
-'0 



V£-.2(l-T2)+r2 + l^ 



jiTM _ gr,2- V^2(l-T2)-Fr^" 

2t 



A-2 



i-TM 



^£^2(1 -r2)+T2' 
gr,2(l"T2) 



2t2 



(er,2(l - t2) + t2)3/2 ^ ^^ ^(l _ r^) + ' 



' Vg.,2(l-T2)+r2(£,,2 + t2) ' 
<2(1--^^)^ ^^(-£r,2^^+gr.2+gr,2 + l) 

^ (£,,2(1 - r2) + r2)3/2(e^ 2 + r^) (srM - r^) + T^){er,2 + t2) 



+ 



(£r,2Vgr,2(l-T2) + r2 + l)' 



(£.,2(1 - r2) + r2) (v/£r,2(l-r2)+T2 + £,,2) ' 
Notice that TJ,/C^,/C2 only depend on £^,2 and r. They are independent of li and e. 
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Gathering the expansions obtained above, we can write 

^{l+li)m,il+li+i)m 

1 /T,™ T fi-r\''+'''-'-^-^ r fm^ 21 16^ 2el „ , 

dOexp — \-{k + k+i)0 

' Z T T T 



.2fe„,_/ ,2 , 2iTO^\ /""^ , , ,2fc „..„ ,_/2 , 2im(^' 



exp (^-^^ + — ^ j y d^V"" exp (^-^'^ + — ^ j (1 + + 0,,2) 

1 + A,2 Ml WTo™(i + 0/C™ + 02^te^ Wl+A/;,2 All 

A^i 1 + A,2 j V To™ (1 + + j 1, A^i 1 + ^,,2; ' 

where 

exp + B,^2 + 2^i,l + 2:'i,2 + ^i,l + ^i,2 + Ui^l + "^^,2) (1 + A,2) (1 + Ci,2) (1 + + <?i,2) (1 + Qi,l + ei,2) 

Notice that A^i,/C|,/C2 ^■'^^ independent of k,k',ip,ip' and ^. Performing the summation over k and k' using the 
formulas 

OO V. 

^ fc! ' 

OO u 

V 



fc=0 

OO ^ 

fe=o 



we obtain an expansion of the form 

M,+,)„,,+,,0™ ( _j.TM j j y_ exp ------ + ft + k^,)Oj 

^Tq™ (1 + + + Te^) + f™Ml (fo™ + To™) 

(fo™ + fo^^i) 1 To™ (1 + ^/C?M + ^2/c™ + Tea) + To™A^?^ 

The TZ2 term comes from £j 2 and A/i,2- The Gaussian integrations over ip and (p' can be performed straightforwardly 
and give 

l/T™. \ 1 + r fl-T 



Hl+U)m,{l+U+,)m -T™J ^ Vl+^. 

X d^exp (-^ - ^ - ^ - ^ + (Z, + ;,+i)0) (1 + Qi,i + Q,,2) 

0™ (1 + + + ^^2) + fo™A^ f (fo™ + To^m) 1 

(fo™ + To™) Ml To™ (1 + eJCj"^ + e^/CjM + U2) + To™At2^ 

ZY2 comes from TZ2 and it is independent of 6. Before performing integration over 6, one is supposed to multiply 
(1 + Qi,i + Qifi) into the matrix after it. Up to the terms of order e, we can write 

1 + Qi,i + Q.i,i « (1 + Qi,2) (1 + Qi,i) , 
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and only multiply (1 + Qi,i) into the matrix. On the other hand, up to the terms of order e, we can extract the term 
U2 of order e out from the matrix. These give 



{l+li)m,(l+li+x)m 

ii+'i+i 

n \ / 1 _ -^N 

l+li 



n T^TM 111, 



J rf^exp ( _±_l^_f£i + (Z. + l,^,)e\ (1 + +Z^2) 



(fo™ + To™) TWi To™ (1 + Qi,i + + 0Qi,i/C™ + 0^IC™) + f^'^Mj^ 

Performing the Gaussian integration over 6, we have 

„..W \ exp + 77(«i + «i+i)' (l+^i,2+W2) 



^fo™ (1 + + V™ + + v™) + fJ^Ml (fo™ + To™) 

(fo™ + To™) All To™ (1 + + + 5i,2/C™ + V^^m) + fJ^Ml^ 



where 



(10) 

n poo / jn2 \ 

Sij J d0exp (-— + ik + k+i)0\ Qi,j, j = 1, 2, 
5,,2 / rf^exp + {k + k+i)0 0Qi^i, 

V-2 = (j + ^a. + /m)^)/C^ 

Next we consider the small e expansions of • Debye asymptotic expansions of modified Bessel functions say 
that: 

r / N 1 e'"'(^^ / ui(r(z)) 



where 



2zy(l + z2)l/4 \^ J, 



t(0) = J-^ , V{z) = Vl + z^ + In 



Let 



UJ 1 

J-, 2;i=ni2:, i^ = l + li + -. 
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Then we find that up to terms of order e, we have 



(zi)) ™i(t(z)) 



In small e expansion, 

(i^) exp - + lU)+X,,, +X,. 

Therefore, we have an expansion of the form 



-T^T^y- I 



/ 1 - T \ — ^ — ^ / 2/ r 

(^Y^j exp - + 2i,2 

To™(l + j;TM + ^TM^j, 



where 



iTE 




V£r,i(l-r2) + r2 + l' 



TM - i/er,l(l - T^) + 







£r,i + v'£r,i(l-r2)+r2' 



with 



4r 



V£r,i(l-T2)+r2- 



•^2 — — ; T- — /o 1" 



,TM 



3^2 



(£,,l(l - r2) + t2)3/2 ^2 _ ^2) + ^2) (1 _ ^2) + ^2 + 

4(l-r2)(r^+c.,i) 
r2 (£,,i(l-r2)+r2)' 

r _ 8£r,ir2 - 3£r,i - Ser.ir'' + Qr^ + Sr"^ 

(e,.i(l-T2)+r2)i/2 12(£^,i(l-T2) + r2) ' 



V£r,l(l-T2)+r2(T2 +£,,!)' 

£,,i(l - r2)(8T2 +4t4 + 4£,,i - 4e,,ir4) 4(1 - r^fel^ (l + ^ SrA^ - t'') + t^)' 



{er,i + r^)(.M(l - r^) + r^f^' ,2 _ ,2) + ,2) (^ ^,,(1 - r^) + + £,,)^ 

4<i(l-r2)^ 



(11) 



(t2+£,,i)(£,,i(1-t2)+t2)' 

TM £r,i(l - T^)r 7£^^ir^ - 4£^^it2 - 3£^^i - 5£^,it6 + 13£^,iT^ - 18£^,ir^ + Sr^ - 3r^ 

' ~(e.,i+T2)(£,,i(l-r2)+T2)i/2 ' ' I2(£,^i + r2)(£,,i(l-r2) + r2) 

Notice that Tg , Wi, W2,!y2 only depend on £^,1 and r. They are independent of and e. 



10 



Substituting dTTI) into (|10|) . we have an expansion of the form: 



2 ~ Ml 

^yTEyTE^TE ^ rTEyTM^2 yTE |^yTE ^ fTM^ j^^ 

J.TM f fTE ^ fTM\ ^ ^ yTM^TM ^ T™fj'' Ml 



where 



7i,i —2^1, 



1 



A* =1 + Jiti + 5m + V*i + J,*iV*i + J,:A,i + Ji:2 + + V*2. 

Substituting (jl2p into ([5]), and extracting terms up to order e, we have 

he 



Cas 



°° 1 -1 ^1 _(s+l)/2 ^oo /.oo / ^ /-oo 

1- dr \ / dll-^^-^y^ dm\T\ dh 



X exp 



m2(s + l) 2e/(s + l) 



1=0 



*=TE,TM 

E h^i'MEE^M^^i + EE^M^.i+E'^.iKi+E^^^^^ 



1 + E E ^^.l^J.l + E ^^'2 + {S + l)U2 I + XMl 

2=0 j=0 i=0 



*=TE,TM 



i=0 j"=0 



i=0 jT^i 



i=0 



i=0 



where 



and 



X =(,s + 1) <^ (To^To™ + To™ro™ 



TnTET^TE 

^0 ^0 



T^TMo^TM 

^0 -'o 



s+1 



j-TEj-TE _ jjTM jjTM 



QrynTErrnTErrnTMrriTM 

^-'o -'o -'o -'o 



t^TEt^TE 
-'O -'o 



(12) 



E E ^ti^i+E E Ki^i+E^:iE^i+E^:2+EK2i 

i=Oj=i+l i=Oj"=i+l i=0 j=0 i=0 j=0 



-Zi2 =7i,2 + 7i i5i 1, 



^^TM^^TM 

-'o -'o 



T^TEtiTE _ rpTMrpTM 
^0 ^0 ^0 ^0 



We have omitted those terms of order -y/e since they are odd in one of the k and thus would give zero after integration 
with respect to li. It follows that 

Eq^^ is the leading order term that comes from those terms of order e", and i?cas next-to- leading order term 

that comes from those terms of order e. e^,! and £^,2 are functions of 



R T 

They are independent of li and m. Performing the Gaussian integration over h, 1 < i < and to, we find that 



_ he ^ 1 



4^^^(s + l)'7o tVT^Jo 



dT 



dl I exp 



2eZ(s + 1) 



(13) 



*=TE,TM 



s=0 V I / ' V J- ' \ / l^*=TE,TM 

The exphcit formulas for , and Si* are given by 

^ "IT ^^"^ ^ "^^^ ^ ^ "^^^ ^ i ^^^^ ~ ^ "^^^ ~ ^^^^ ~ "^^ ' 
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(14) 



YIIt ^ • ' • 21t 31 s + r 

l-r2 



2/r(s + 1)' 

=Dyy/Cf + DvjIClW*, + DjjWf + +Dv]lCl+ DjW* + (s + 1)^3^* 



with 



Cv 


= -f((. + l)-^ 


+ 2(. + l)) + i^(,s + l)2 + 


2/(^ + ^) + 


Cj 


= -f((« + l)^ 




5 


Dvv 


= l5l((^ + ^)^- 


2(s + l)2 + 2(s + l)-l), 




Djj 


= 4^1 ((^ + ^)^ - 


2(s + l)2-(s + l) + 2), 




Dvj 


-1^((^ + ^)^- 






Dv 




3(s + l) + l), 




Dj 


= lil((^ + ^)^- 







121 



Using the fact that Tq*,To*,/Ci,/C2, ^i, ^2,3^2 are independent of e, it is straightforward to take derivative with 
respect to d. For the Casimir force 

p _ dEcas 

we find that 

where Fq,^ and i^^^^g are respectively the leading order and next-to-leading order terms with 
2nR^^s + lJo T^VT^Jo '^"^''PV r ) ^ 



12 



Here 



(G"^ + 1)' + 2Gs + 1)) - I (2(s + If + 3t(s + 1) + 1) , 

-r'' + 5t2 - 12 , , 1 + r - r2 t 1 
(s + 1)- 









er 




~ ~ y 




er 





12Zr V ' ^ ■ 2/t 6/s + l' 

.s + l)^ + 2Gs + l)) + l(,s + l)^ + J(. + l) + 2_, 



For the force gradient dFcas/dd, the leading order and next-to- leading order terms are 



s+1 



where 



((s + 1)3 + 2(,s + 1)) - ^ ((2 + r2)(s + 1)2 + 3r(s + 1) + 1 + 2r2) 



(.s + 1) + 



12;t ' ' 21t 
^=-?((^ + lf + 2(. + l)) + i^(.s + l)^ + ^(. + l)+' + ^"' 



3 ' / ■ \ ■ // ■ 6/ ' ' 2/' ^ 12Z 
5} = -^((. + 1)3- (.s + 1)) + ((. + 1)^-1). 



Let us compare the leading order term to the proximity force approximation. The Casimir energy density between 
a pair of parallel dielectric plates with relative permittivities er,i and er,2 is given by the Lifshitz's formula [2ll |: 

t /"OO /"OO 

£Lid)~J dK dqq ln(l-rJr*e-2'.''), (15) 

Jo Jk *=te,tm 



where 



y^{er,i - 1)^2 + (f -q 
^{er,i - 1)k2 + g2 + g,' 
^TM _^r,i1 ~ \/ i^r,i ~ 1)^^ + (/^ 

The proximity force approximation to the Casimir interaction energy between a sphere and a plate with relative 
permittivities e^,! and £^,2 is given by 

/•ao 

E^l^ =2nR / du4^^{u). (16) 

Jd 

Expanding the logarithm in (fT5|) and substitute into (IT6)) . we find that 

, p CXj -| />00 /"OO /"OO 



s=0 \ ■ / .'K, *=TE,TM 
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Now making a change of variables 



I iVl-T 



2 



1 = TT' 



Rt' Rt ' 



we finally obtain 



where 



„TE 



V£^,,(l-r2) + r2-l 
"V£-,,(l-r2)+T2 + l' 

Compare to (I13p . we find that our result for the leading order term agrees completely with the proximity force 
approximation. 

IV. PLASMA MODEL 

In this section, we consider the special case where the dielectric permittivities of the sphere and the plate are 
described by the plasma model: 

= 1 + ^, 

where ujp^i is the plasma frequency of the material. 
Let 

c 

In terms of the variables 

el 

t = — 

T 



and T, we have 

First consider the case 

Wd,i > 1. 

The limit where uJd,i — ^ oo,i — 1,2, is the perfect conductor limit. We can compute analytically the asymptotic 
expansion of the leading and next-to-leading order terms in the small parameters 

1 

tti = . 



Specifically, we have 



3^ ^ / oo oo ^ oo oo 

Ec.s « - E E + ]^ E E ^^^«>2 + ■••], (18) 

\ i=0 j=0 i=0 j=0 
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TABLE I: The coefficients 



p 


exact value 


numerical value 


A),o 


1 


1 


/3i,o 


4 
~3 


— 1.3333 




9 

5 


1.8 




18 
5 


3.6 


P3,0 


16 32 2 
y + 735" 


— 1.8560 




48 
7 


-6.8571 


/34,0 


25 326 2 
9 1323" 


0.3458 


P3,l 


100 326 2 
9 1323" 


8.6791 


flo o 


50 
3 


1 6 6667 


/3s,o 


36 1220 2 379 4 

\ TT TT 

11 1617 32340 


3.0322 


/34,1 


180 2440 2 

1 TT 

11 1617 


-1.4707 


/33,2 


360 1220 2 
11 + I6I7" 


-25.2808 



where /3o,o = 1 a-nd 

pPFA.PC ^ TT^hcR 

is the leading order approximation to the Casimir interaction energy between a perfectly conducting sphere and a 
perfectly conducting plate. The exact values of f5i,j and for i + j < 5 are listed in Table U and Table HIl From 
(|13p. it is obvious that the leading term is symmetric when we interchange £r,i with £^,2- It follows that 

Pi,j=P],i foraU(i,j). 

Hence, we only list the coefficients of /3ij when J > j in Table ID 
From (fT8|) . we have 

E E + I E E A.„ala^ + . . .y (19) 



dF, 



Cas 



dd 



TT^hcR 
120# 



CX3 00 

EE 

=0 j=0 



(i + j + 2)(z+j+3) 
6 



fli fto + — 



^>F^(*+J + l)(* + J + 2) 



* J 1 2 



(20) 
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TABLE II: The coefficients Ai. 



exact value 



numerical value 



Ao,o 
Ai,o 
Ao,i 

A2,0 

Ai,i 

Ao,2 

As.o 

A2,l 
Al,2 
Ao,3 
A4,0 
As,! 
A2,2 
Al,3 
Ao,4 

As.o 

A4,l 
A3, 2 
A2,3 
Al,4 
Ao,5 



_20 1 
"^"^ 3 

56_1_ 
3 7r2 

56_I_ 

3 7r2 

398 1 
" 21 7r2 



32 
45 

14 
45 

401 
315 

454 



796 1 

" 21 7r2 315 



398 1 113 
"ir^ ~^ 315 



410 1 
'21 



37 286 2 

1 TT 

18 6615 



410 1 
7 7r2 

410 1 



26 
7 

16 
7 



410 1 

"2r^ 



126 6615'^ 



69824 1 35141 
+ 



3465 7r2 

279296 1 
3465 7r2 

139648 1 



+ 



10395 

84176 
10395 



28022 2 
99225 ^ 

2774 2 



14175 

32 . 



742 

1155 7r2 "99"^ 11025 



279296 1 43856 
+ 



3465 7r2 
69824 1 



+ ■ 



10395 
14981 



46558 2 
1091475^ 



11962 2 

TT 



3465 7r2 ■ 10395 1091475 

26732 1 150368 



1287 7r2 



4937399 2 

-TT ■ 



27027 5675670 



1142 4 

TT 



63063 



133660 
1287 


1 

7r2 


35026 773884 
2079 ' 567567 


267320 
1287 


1 

7r2 


548024 
27027 


26212 
51597 


267320 
1287 


1 

7r2 


415724 
27027 


16826 
81081 


133660 
1287 


1 

7r2 


256888 
27027 


19984 
81081 


1 84218 


3329 


2 , S 



1287 7r2 27027 62370 



8059 4 
7r 



2522520 



-1.6931 



1.1802 



1.5802 



-0.6473 



-2.3993 



-1.5615 



0.3493 



2.2202 



3.6488 



1.3527 



-1.4484 



-2.0007 



-4.7269 



-4.3690 



-0.7087 



3.3627 



7.1324 



5.7822 



7.7116 



3.4503 



-0.1736 
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Here 



and 



PFA.PC 



Cas 



dF, 



PFA.PC 

Cas 



dd 



TT^hcR 
360^3 



TT^hcR 
120^4 



are respectively the leading order approximations to the Casimir force and force gradient between a perfectly con- 
ducting sphere and a perfectly conducting plate. Setting ai = 02 = in , and (PU)) , we obtain 



Ecas 

Fcas 
dFcas 

dd 



iE, 



PFA.PC 
Cas 

PFA,PC 
Cas 

PFA.PC 



1 + 



1 + - 



dF, 



Cas 



dd 



1 d 



20" 






+ 


10' 






+ 


20 " 


3 


7r2 



which are well-known results for the leading and next-to-leading order terms of the perfectly conducting sphere-plate 
configuration [HI, [T^. 

Next we consider numerical results with uip^i — Wp_2 — 9eV, which is the plasma frequency for gold '22']. The radius 
of the sphere R is taken to be 1mm. Substituting (jl7p into the formulas obtained in the previous section, we can 
compute numerically the leading order term (the proximity force approximation) and the next-to-leading order term 
of the Casimir interaction. In Figs. [TJ [2]and[3l we plot the leading order term, the sum of the leading order and next- 
to-leading order terms of the Casimir interaction energy, Casimir force and force gradient, normalized respectively by 
-^Casf ^^' -^Cas^'^^ ^'^'^ ^^cls^'^'" / ^ fuuctiou of d for d between lOnm and 100/xm. From the figures, we notice 
that when d/R ^0.1, the corrections to PFA become significant and they would contribute corrections of about 10%. 




FIG. 1: The leading order term of the Casimir interaction energy normalized by S^fif '^'^ (dashed line) and the sum of the 
leading and next-to- leading order terms normalized by E^^^'^''^ (solid line). Inset is the ratio of the latter to the former. 



To have a better picture about the corrections to the proximity force approximations, define 0i,_b, Oi^f and 9i by 



yi.E 



R^Cas 

'dF^^J 



^R dF^Jdd 
' ddF^Jdd' 
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FIG. 2: The leading order term of the Casimir force normahzed by Fq^^'^'^ (dashed hne) and the sum of the leading and 

PFA PC / 

next-to-leading order terms normalized by Fq^^ ' (solid line). Inset is the ratio of the latter to the former. 




FIG. 3: The leading order term of the force gradient normalized by OFq^^'^''' /dd (dashed line) and the sum of the leading 



and next-to-leading order terms normalized by dF^ 



/dd (solid line). Inset is the ratio of the latter to the former. 



so that dJ) and ^ hold. 

In Figs, m we plot 9i^e, 6*1, f and 9i as functions of d for d between lOnni and lOO/xm. As d increases to 100/im, 
we find that 6'i,e, 9i^f and 9i tend respectively to the values 1/3 - 20/71^ = -1.6931, 1/6 - IO/tt^ = -0.8465 
and 1/9 — 20/(37r^) = —0.5644, which are corresponding values for perfect conductors. This is not surprising since 
uj(i,i — ojp^id/c is approximately equal to 5000 when d = 10~*. When d is small, the deviations from the limiting 
values for perfect conductors are very significant. On the other hand, we also notice that 9i^f and 9i are bounded 
below. 9i is a quantity that can be measured experimentally [5|. From Fig. 21 we find that it is bounded below by 
-0.57. 



V. DRUDE MODEL 



The Drude dielectric function is given by 



w2 . 



where 7i is the relaxation frequency of the material. In the limit where 7i — >■ 0, the Drude dielectric function becomes 
the plasma dielectric function. 
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Let 



In terms of the variables t and t, we have 



lid 



ld,i = • 

c 



Substituting this into the results obtained in Section [llll we can compute numerically the leading order and next-to- 
leading order terms of the Casimir interaction for Drude models. 

Let us consider the case where Wp^i = 2 = 9eV and 71 = 72 = 0.035 eV, which are the conventional values used 
for gold [22} . 





E°, Dmde 






1 


e", plasma 






E°+E', Drude 






0.8 


^ E^+E', plasma ^ 






0.6 


J' 

/ 0.95 






0.4 


/ 0.9 

/ 0.85 


Drude \ 

<> plasma \ 




0.2 








10 


, 10:' 10^ . 



10"** 10"' ,Jof 10^ 10"* 
a{ra) 

FIG. 5: The leading order term and the sum of the leading order and next-to-leading order terms of the Casimir interaction 

PFA. PC 

energy normalized by ' . The inset shows the ratio of the latter to the former. 



In Figs. [5l[6]and[71 we plot the leading order term, the sum of the leading order and next-to- leading order terms of 
the Casimir interaction energy, Casimir force and force gradient, normalized respectively by E^^^'^'^ , Fq^^'^'^ and 
dF^^^'^"^ / dd, as a function of d for d between lOnm and 100/im. Both the Drude model and the plasma model are 
plotted on the same graph to show the comparison. To get a better picture, we plot the ratio of the plasma model 
to the Drude model for the sum of the leading order and next-to-leading order terms in Fig. [S] From the figure, we 
notice that if the plasma model is used instead of the Drude model, the error is at most 2%. 

In Figs, m we plot 9i^E,di.F and 9i for Drude model and compare to that for plasma model. As for the plasma 
model, we notice that for the Drude model, as d increases, 9i^e, 9i,f and 9i tend respectively to the limiting values 
1/3 - = -1.6931, 1/6 - IO/tt^ = -0.8465 and 1/9 - 20/(37r2') = -0.5644, the corresponding values for perfect 

conductors. When d is small, the deviations from these limiting values are very significant. On the other hand, 9i,p 
and 9i are also bounded from below. 
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F°, Dmde 




FIG. 6: The leading order term and the sum of the leading order and next-to-leading order terms of the Casimir force 

PFA. PC 

normalized by F^^^ ' . The inset shows the ratio of the latter to the former. 




FIG. 7: The leading order term and the sum of the leading order and next-to-leading order terms of the force gradient 

PFA. PC o 

normalized by a-Fc^s ' 9d. The inset shows the ratio of the latter to the former. 

The ratios of the plasma model to the Drude model for Oi^e,(^i.f and 6i are plotted in Fig. (TU] From the figm'e, 
we find that if the plasma model is used instead of the plasma model, the error is not more than 4.5%. 

VI. CONCLUSION 

Starting from the functional determinant representation of the Casimir interaction energy, we have used the per- 
turbation method to obtain analytically the leading order and next-to-leading order terms of the Casimir interaction 
energy, Casimir force and force gradient for the interaction between a sphere and a plate. The results are written as 
double integrals over functions of the dielectric permittivities of the objects, and are hence general. The leading order 
terms are shown to equal to that predicted by proximity force approximation. The results on the next-to-leading 
order terms are new, and they complement those obtained in [l^ using derivative expansion. 

With given dielectric permittivities of the sphere and the plate, the double integrals representing the leading order 
and next-to-leading order terms can be computed numerically, and this is demonstrated for a gold sphere in front of 
a gold plate, where both plasma and Drude models are used for the dielectric functions of gold. It is observed that 
even at d/R ^ 0.1, the next-to-leading order term would contribute a correction to the leading order term of about 
10%. Of particular interest is the ratio of the next-to-leading order term divided by d/R to the leading order term, 
denoted by 9i. It is found that when cupd/c is large enough, 9i{d) tends to the corresponding limiting value for perfect 
conductors. However, when Wpd/c is small, the deviation from the limiting perfect conductor value is significant. This 
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1.02 




FIG. 8: The ratio of the plasma model to the Drude model for the sum of the leading order and next-to-leading order terms. 




10* 10 ' ,JOt 10 ' 10"* 10 ' 10 ' ,10? 10 ' 10"* 10* 10 ' ,10? 10 ' 10"* 

d (m) d (m) d (m) 

FIG. 9: 6\^E{d), 9i.F{d) and 6i{d) computed using Drude model (solid line), compared to those computed using plasma 

model (dashed line). 



signifies that in the nano range, we cannot model real metals by perfect conductors. 

A comparison between plasma model and Drude model shows that their difference is below 2% for the sum of the 
first two leading order terms and below 4.5% for the values of 6i{d). In fact, this small difference is expected at zero 
temperature. In this work, we haven't considered the thermal effect. When d is small enough such that 27rkBTd/hc 
is ^ 1, thermal effect can be neglected. For example, when T = 300K, thermal effect can be neglected when d <C 
1/im. Nevertheless, it would be interesting to examine the behavior of the Casimir interaction at room temperature 



1.05 




10-^ 10-' ^(1^-; 10- 10-^ 

FIG. 10: The ratio of the plasma model to the Drude model for 6i,e, di,F and 8i. 
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and its interplay with material properties. This would be addressed in a forthcoming work. 
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